Critical exponents of the quantum phase transition in a planar antiferromagnet 
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We have performed a large scale quantum Monte Carlo 
study of the quantum phase transition in a planar spin-1/2 
Heisenberg antiferromagnet with CaV409 structure. We ob- 
tain a dynamical exponent z = 1.018 ± 0.02. The critical ex- 
ponents /3, V and -q agree within our errors with the classical 
3D 0(3) exponents, expected from a mapping to the nonlin- 
ear sigma model. This confirms the conjecture of Chubukov, 
Sachdev and Ye [Phys. Rev. B 49, 11919 (1994)] that the 
Berry phase terms in the planar Heisenberg antiferromagnet 
are dangerously irrelevant. 

Instead of classical transitions controlled by tempera- 
ture T a quantum phase transition between a symmetry 
broken phase with long-range Neel order and a quantum 
disordered state with a finite spin excitation gap may be 
realized at T = by controlling a parameter g to in- 
crease quantum fluctuations. Criticalities around such 
quantum phase transitions a,t g = g^ may reflect inher- 
ent quantum dynamics of the system and yield unusual 
universality classes with rich physical phenomena. 

The most prominent example are the high temperature 
superconductors. There the quantum spin fluctuations 
are thought to lead to d-wave superconductivity as soon 
as antiferromagnetism is suppressed by hole doping. 

In this letter we want to discuss a two dimensional 
quantum Heisenberg antiferromagnet (2D QAFM) that 
exhibits such a quantum phase transition. We will 
present strong numerical support for the conjecture that 
the 2DQAFM is in the same universality class as the 
quantum nonlinear sigma model (QNLcrM) even when 
Berry phase terms are present. 

The universality class is characterized by the critical 
exponents. Approaching the quantum critical point from 
the disordered side the spatial correlation length diverges 
with the correlation length exponent v. The space and 
time dimensions are however not necessarily equivalent, 
and the correlation length in the time direction diverges 
in general with a different exponent zv, where z is the 
dynamical exponent. In a Lorentz invariant system space 
and time directions are equivalent and z = 1. Related to 
the divergence of the correlation length is a vanishing 
of the spin excitation gap with the same exponent ziy. 
When passing through the critical point long range order 
is established. The order parameter in the case of a Neel 
ordered antiferromagnet is the staggered magnetization 
nis- Near the critical point tti^ vanishes with the order 
parameter exponent /3. At the critical point itself the 
real space staggered spin correlation shows a power-law 
falloff with power 2 — d — z — rj, where rj is the correlation 
exponent. These three exponents are related by the usual 



scaling law 



2(3={d + z-2 + Tfjv, 



(1) 



z in a quantum sys- 



where the effective dimension is d - 
tern. 

Quantum critical behavior of a planar antiferromag- 
net has been intensively studied by a number of groups. 
Most analytic calculations are based on the QNLcrM. 
The critical exponents of the QNLcrM can be determined 
from simple symmetry, universality and scaling argu- 
ments |l|,yi. Lorentz invariance implies that z = 1. Fur- 
thermore the 2D QNLcrM is equivalent to the 3D classical 
sigma model. This in turn is in the universality class of 
the 3D classical 0(3) model, or the classical 3D Heisen- 
berg ferromagnet. The exponents /3, v and rj should thus 
be the same as the well known classical exponents of these 
models (see Tab. I). 

Chakravarty, Halperin and Nelson have discussed the 
phase diagram of a planar Heisenberg antiferromagnet by 
using the QNLcrM. In their discussions they concentrate 
on the ordered phase and describe it as a classical 2D 
antiferromagnet with renormalized parameters. 

Chubukov, Sachdev and Ye have investigated the quan- 
tum critical regime of the QNLcrM in close detail. They 
make some further predictions based on scaling argu- 
ments. On the ordered side the spin stiffness ps vanishes 
as 



Ps oc {gc - g) 



{d+z-2)u _ 



{9c - gY, 



(2) 



where the second equivalence comes from the prediction 
that z = 1. They also predict that the uniform suscepti- 
bility at the critical point is universal: 
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Here c is the spin wave velocity and rii(oo) a universal 
constant. Estimates for Sli(oo) are listed in Tab. II. 
The spin wave velocity c scales as 



c^ig^g^r^^-'^ 



(4) 



and is thus regular at the critical point if z = 1. 

The equivalence of the 2D QAFM to the 2D QNLcrM 
however is still an open question because of the existence 
of Berry phase terms in the QAFM that are not present 
in the QNLcrM M. It has been argued that these terms 
cancel in special cases, such as in the bilayer model [Q,||. 
Then it is plausible that the quantum phase transition 
is in the same universality class as the QNLcrM. This 



was confirmed by quantum Monte Carlo calculations of 
Sandvik and coworkers |^,0,||- They have investigated 
the finite size scaling of the ground state structure factor 
and susceptibilities on lattices with up to 10 x 10 x 2 
spins. Although these lattices are quite small they still 
found good agreement of the exponents z and t] with the 
QNLcrM predictions |^,0 (see Tab. I). In another study 
Sandvik et al. M have investigated finite temperature 
properties of the bilayer QAFM on larger lattices and 
also found good agreement with the QNLcrM predictions. 
In the absence of Berry phase terms the equivalence of 
the QAFM and the QNLcrM is quite well established by 
these simulations. 

But in general these Berry phase terms exist. 
Chakravarty et al. argue that they can change the 
critical behavior and lead to different exponents |l|,[lO|. 
Chubukov et al. on the other hand argue that the Berry 
phase terms are dangerously irrelevant and do not 
influence the critical behavior. Previous numerical simu- 
lations on dimerized square lattices |^,|ll| are indeed not 
consistent with the QNLcrM predictions. The reliabil- 
ity of the results however is questionable because of the 
restriction to very small lattices of 12 x 12 spins and be- 
cause of complications with scaling arising from inequiv- 
alent spatial directions. On the other hand the discrep- 
ancy could be an effect of the Berry phase terms that are 
present in the dimerized square lattice but probably not 
in the bilayer. 

Using the new quantum cluster algorithms 1 12 13] we 



could simulate much larger lattices at lower tempera- 
tures. On these larger lattices we find perfect agreement 
with predictions made based upon the 2D QNLcrM de- 
spite the presence of Berry phase terms. 

As the universality class of a phase transition does not 
depend on the microscopic details of the lattice structure 
we are free to choose the best lattice for our purposes. 
We have chosen the CaV409 lattice, a 1/5-th depleted 
square lattice depicted in Fig. 1 for our calculations. 
There are three reasons for this choice. Firstly the Berry 
phase terms are present on this lattice [Q. Next both 
space directions are equivalent, in contrast to the dimer- 
ized square lattice ||^,|ll| . This makes the scaling analysis 
easier. Finally at the quantum critical point all the cou- 
plings are nearly equal in magnitude, which is also opti- 
mal from a numerical point of view. We have performed 
our simulations on square lattices with N = Sn^ spins, 
where n is an integer. Our largest lattices contained 20 
000 spins. For the following discussion it is useful to 
introduce the linear system size L in units of the bond 
lengths a of the original square lattice: L = y/5N/Aa. 

The phase diagram of this lattice has been discussed 
in detail in Ref. |15|. By removing every fifth spin we 
obtain a lattice consisting of four-spin plaquettes linked 
by dimer bonds. We label the couplings in a plaquette 
Jo and the inter-plaquette couplings Ji. By controlling 
the ratio of these couplings Ji/Jq we can tune from Neel 



order at Ji = Jq to a quantum disordered "plaquette 
RVB" ground state with a spin gap A = Jq at Ji = 0. 

At some intermediate coupling ratio {Ji/Jo)c the sys- 
tems has a quantum phase transition. The first step in 
the determination of the critical behavior is a high pre- 
cision estimate of the critical point gc- We have calcu- 
lated the second moment correlation length ^l on sys- 
tems of various sizes L. The temperature was chosen 
to be ksT = Jga/ L, keeping the finite 2 -f- 1 dimen- 
sional system in the cubic regime. From standard fi- 
nite size scaling arguments it follows that this correla- 
tion length ^L scales proportional to the system size L at 
criticality. We have calculated the ratio £,l/L (shown in 
Fig. 2) for a variety of couplings and system sizes up to 
N = 9600 and have determined the critical coupling to 
be (Ji/Jo)c = 0.939 ±0.001. 

Next we have calculated the finite size scaling of both 
the staggered structure factor S'(Q) = L^nis and of 
the corresponding staggered susceptibility. At critical- 
ity they scale like 
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(5) 
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The temperature was chosen to be ksT — Joa/(4_L). 
This was low enough to see the ground state proper- 
ties on the finite lattice. By fitting our results shown 
in Fig. 3 we obtain the estimates z — 1.018 ± 0.02 and 
77 — 0.015 ± 0.020. This is perfectly consistent with the 
Lorentz invariance (z = 1) expected from a mapping to 
the QNLcrM. We will discuss rj below together with the 
other exponents. From these fits it is also obvious that at 
least N ~ 800 spins are necessary to obtain good scaling. 
The remaining exponents /3 and v are best calculated 
from the magnetization rris and the spin stiffness ps on 
the ordered side. Good estimates for rris and ps can 
be obtained from the Hasenfratz-Niedermayer equations 
||l6| . These authors have calculated the exact finite- 
size and finite-temperature values of the low-temperature 
uniform and staggered susceptibilities Xu and Xs for the 
ordered phase of a 2D QAFM on a lattice with the sym- 
metries of a square lattice. Their equations, determined 
by chiral perturbation theory, are correct for the low 
temperature regime ksT <^ 'ii^Ps with cubic geometry 
ksTL/fic K, 1. Up to second order in T (or 1/L re- 
spectively) the susceptibilities are universal, determined 
by only three parameters: the staggered magnetization 
TOs, the spin stiffness ps and the spin wave velocity c. 
Two high precision auantum Monte Carlo studies have 
confirmed their equations for the square lattice QAFM 

We have calculated the susceptibilities for a wide range 
of couplings 0.95 < Ji/Jq < 1.1, lattice sizes 800 < N < 
16200 and temperatures 0.006 < T/Jq < 0.1. The fits 
to the Hasenfratz-Niedermayer equations are all excel- 
lent, with x^/d.o.f. « 1.5. This is another confirmation 



of the universality of the Hasenfratz-Niedermayer equa- 
tions. From the fits we obtain the staggered magnetiza- 
tion rus, the spin stiffness ps and the spin wave velocity 
c. The exponents (3 and v can then be obtained in a 
straightforward way (see Fig. 4) and are listed in Tab. I. 

Let us now discuss the results. First we observe that 
the exponents satisfy the scaling relation Eq. ffl), con- 
firming the validity of the scaling ansatz for this quantum 
phase transition. The exponents /3, v and rj are in excel- 
lent agreement with the exponents of the 3D classical 
0(3) or Heisenberg model. They are however incompat- 
ible with the mean field exponents suggested by Katoh 
and Imada from their calculations on small lattices. 

Assuming Lorentz invariance {z = 1) we can improve 
our estimates for the other exponents. The agreement 
of the improved estimates with the 3D 0(3) exponents 
becomes even better. We can rule out not only the mean 
field universality class suggested by Ref . [l|] , but also the 
Ising universality class (see Table I for a comparison). 

This excellent agreement is a strong numerical sup- 
port for the conjecture of Chubukov, Sachdev an Ye ||] 
thet the Berry phase terms in the 2D QAFM are indeed 
dangerously irrelevant. To further confirm their predic- 
tions we have calculated the uniform susceptibility close 
to criticality down to T = 0.02, more than an order of 
magnitude lower than Ref. . We have extrapolated the 
finite size results on lattices with up to iV = 20000 spins 
to the thermodynamic limit. Looking for the coupling at 
which a linear behavior occurs gives an independent es- 
timate of the critical point: ( Ji/ Jo)c = 0.939 ± 0.002, in 
excellent agreement with the above estimate. The linear 
slope is ni{oo) {Jo /hc)'^ = 0.238 ± 0.003. By extrapo- 
lating the spin wave velocity determined in the ordered 
phase by the Hasenfratz-Niedermayer fit to the critical 
point we get he/ J = 1.04 ± 0.02 and thus rii(cx)) = 
0.26 ±0.01, again in excellent agreement with Chubukov 
et al. H (see Tab. II). 

To summarize, we have calculated the critical expo- 
nents of the quantum critical point in a planar antifer- 
romagnet by a large scale quantum Monte Carlo study. 
Our exponents agree perfectly with predictions made by 
a mapping to the 2D quantum nonlinear sigma model. 
The dynamical exponent is z = 1.018 ± 0.02, consistent 
with Lorentz invariance. The other exponents agree with 
the 3D classical 0(3) exponents. The conjecture that all 
quantum phase transitions between a Neel ordered and 
a quantum disordered state in 2D Heisenberg antiferro- 
magnets, whether they contain Berry phase terms or not, 
are in the same universality class as the quantum nonlin- 
ear sigma model is strongly supported. 
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FIG. 1. Lattice structure of the 1/5-th depleted square lat- 
tice of CaV4 09. The dashed square indicates the eight spin 
unit cell used in our calculations. 



FIG. 3. Finite size scaling of the staggered structure factor 
and susceptibility at the critical point. The straight lines are 
fits to the finite size scaling forms Eqns. (pi) and IM). 
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FIG. 2. Plot of the ratio of correlation length divided by 
system size £,l/L. At the critical point the correlation length 
calculated in a finite system is proportional to the system size. 
This is the case for g = Ji/Jq = 0.939 ± 0.001. 



FIG. 4. Staggered magnetization ttis and spin stiffness ps 
calculated by a fit of the low temperature susceptibilities on 
finite lattices to the Hasenfratz-Niedermayer equations [M. 
The straight lines are fits used to obtain the exponents /3 and 



TABLE I. Critical exponents /?, u, rj and z. Listed are both the estimates without making any assumption for z, and the 
best estimate if Lorentz invariance {z — 1) is assumed. For comparison the exponents of the 3D classical Heisenberg (0(3)) 
model, the 3D Ising model and the 2D quantum mean field exponents are listed. The errors given include the uncertainties in 
the critical point. 



model 



P 



V 



2D QAFM 

Lorentz invariant 2D QAFM 



bilayer QAFM || 



0.685 ± 0.035 
0.695 ± 0.030 



0.345 ± 0.025 
0.345 ± 0.025 



0.015 ±0.020 
0.033 ± 0.005 



1.018 ±0.02 
1 (assumption) 



0.03 ±0.01 



1.08 ±0.05 



3D 0(3) 
3D Ising [| 
mean field 



0.7048 ± 0.0030 

0.6294 ± 0.0002 

1 



0.3639 ± 0.0035 
0.326 ± 0.004 

1/2 



0.034 ± 0.005 

0.0327 ± 0.003 





TABLE IL Universal prefactor f2i(cx3) in the linear tem- 
perature dependence of the uniform susceptibility at critical- 
ity. Listed are the results for the quantum nonlinear sigma 
model in a 1/A'^ expansion, the results by classical Monte 
Carlo simulation on a 3D classical rotor model and the re- 
sult of the present study. 



method 



Ref. 



Qi{oo) 



1/N expansion 
classical Monte Carlo 
quantum Monte Carlo 



this study 



0.2718 
0.25 ±0.04 
0.26 ±0.01 



